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Abstract. We report the results of an inviscid spatial stability calculation for a compressible reacting 
mixing layer. The limit of infinite activation energy is taken and the diffusion flame is approximated by a 
flame sheet. Results are reported for the phase speeds of the neutral waves and maximum growth rates of 
the unstable waves as a function of the parameters of the problem: the ratio of the temperature of the sta- 
tionary stream to that of the moving stream, the Mach number of the moving stream, the heat release per 
unit mass fraction of the reactant, the equivalence ratio of the reaction, and the frequency of the disturbance. 
These results are compared to the phase speeds and growth rates of the corresponding nonreacting mixing 
layer. We show that the addition of combustion has important, and complex, effects on the flow stability. 
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1. Introduction. Quite recently it has been realized that an understanding of the stability characteris- 
tics of compressible mixing layers is extremely important in view of the projected use of the scramjet 
engine for the propulsion of hypersonic aircraft For example, Drummond and Mukunda (1988) state: "Even 
though the combustor flow field is quite complex, it can be realistically viewed as a collection of spatially 
developing and reacting supersonic mixing layers that are initially discrete, but that ultimately merge into 
larger more complex zones. These mixing layers begin downstream of a set of fuel injectors that introduce 
gaseous hydrogen in both a parallel and tranverse direction into a supersonic airstream entering from the 
engine inlet. The behavior of the initial portion of the combustor flow, in the mixing layers near the fuel 
injectors, appears to be most critical, since this is where the mechanism for efficient high speed mixing must 
be established to achieve the required degree of combustion downstream. Because of the structure of the 
flow field in this initial portion of the combustor, a single supersonic, spatially developing and reacting mix- 
ing layer serves as an excellent physical model for the overall flow field." Thus knowledge of the stability 
characteristics may allow one, in principle, to control the downstream evolution of such flows in the 
combustor. This is particularly important because of the observed increase in the flow stability at high 
Mach numbers (Brown and Roshko, 1974; Chinzei, Masuya, Komuro, Murakami, and Kudou, 1986; and 
Papamoschou and Roshko, 1986). Because of the increase in stability, natural transition may occur at 
downstream distances which are larger than practical combustor lengths. A number of techniques which 
may enhance mixing are discussed by Kumar, Bushnell and Hussaini (1987). A detailed understanding of 
the linear stability characteristics of compressible reacting mixing layers will be of aid in mixing enhance- 
ment. 

Despite the fact that understanding of the flow field in a reacting compressible mixing layer in a 
scramjet engine is extremely important, there appears to be very few studies of the stability of such flows. 
Menon, Anderson and Pai (1984) studied the inviscid spatial stability of a compressible wake in which there 
was a H 2 -O 2 reaction. These calculations were carried out with a free stream Mach number of 2 and tem- 
perature of 150CFK. When the reaction was turned cm the flow became completely unstable. The phase 
speed was found to be a monotonically increasing function of frequency. It seems that their results show a 
complete absence of neutral or stable disturbances. 

The above result appears to be in conflict with that of Drummond and Mukunda (1988). They carried 
out a numerical simulation using the two dimensional, compressible, time dependent Navier-Stokes equa- 
tions with combustion in a mixing layer. The reaction was the burning of a 10% H 2 , 90% N 2 fuel in air. 
The free stream Mach number was taken to be 2, the temperature above and below the plate was 2000 C, K, 
and the velocities were 2672 m/s and 1729 m/s above and below the plate, respectively. Because of the 
expense and difficulty of carrying out these simulations only a few have been done. Drummond and 
Mukunda found that the nonreacting flow was very stable and that turning on the combustion had little 
effect. But it should be noted that the authors were not carring out a stability calculation, per se, and did 
not excite the flow with disturbances with a fixed frequency. They relied on the "natural" disturbances to 
perturb the flow. 

We have begun a systematic study of the stability of compressible mixing layers in which a diffusion 
flame is embedded. The basic steady flow with which we began is that calculated by Jackson and Hussaini 
(1988). In their study the limit of infinite activation energy was used and the diffusion flame reduced to a 
flame sheet. The flame sheet model is a standard approximation and has been used in the study of the burn- 
ing of a fuel particle in an oxidizing atmosphere and of the flame at the mouth of a tube, for example 
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(Buckmaster and Ludford, 1982; and Williams, 1985). In order to understand the effect of the chemical heat 
release on the stability of this flow, one must first understand the stability characteristics of the nonreacting 
flow. 


In the first part of our study (Jackson and Grosch, 1989, hereafter referred to as Part I), we considered 
the inviscid spatial stability problem for the compressible mixing layer with the mean velocity profile 
approximated by the hyperbolic tangent. We found that there is only a single subsonic neutral mode for 
two dimensional waves, but that there can be three for three dimensional waves. Beyond a critical Mach 
number M„ the Mach number at which the phase speed equals that of a sonic wave, the subsonic neutral 
modes are transformed into supersonic neutral modes which are subsonic at one boundary and supersonic at 
the other (we have not found any neutral or unstable modes which are supersonic at both boundaries). In 
addition, another supersonic neutral mode appears at M . > M„ the Mach number at which the sonic speeds 
of the stationary and moving streams are equal. This supersonic neutral mode has the opposite behavior 
than the previous at the boundaries. That is, if the continuation of the subsonic neutral mode is supersonic 
in the moving stream and subsonic in the stationary stream, this new mode is subsonic in the moving stream 
and supersonic in the stationary stream. Thus, there are always at least two bands of unstable frequencies 
for Mach numbers greater than M*. One of these bands is a group of fast and the other a group of slow 
unstable supersonic modes. The fast modes are supersonic with respect to the stationary stream and the 
slow modes are supersonic with respect to the moving stream. It is important to note that both the fast and 
slow supersonic modes are vorticity modes and neither of them is an acoustic mode (Mack, 1989). These 
groups of unstable modes lie in the frequency bands between zero, corresponding to the sonic mode, and the 
frequency of the supersonic neutral mode. Because these frequency bands always overlap for some range of 
frequencies, there exist two unstable modes at a fixed Mach number and p r (the ratio of the temperature in 
the stationary stream to that of the moving stream) for every frequency in this range. The phase speeds of 
both the fast and slow supersonic modes have a small range about the average, so that little dispersion of 
wave packets is expected, with a reduction in the dispersion as the Mach number is increased. Three 
dimensional disturbances show the same general characteristics as two dimensional disturbances. There is 
always a range of propagation angles for which both the fast and slow unstable modes exist. We also find, 
in agreement with previous studies, that the maximum growth rate for any p r and M occurs for three dimen- 
sional waves. A decrease in p r results in an increase in the growth rate of the unstable waves at any Mach 
number. An increase in the Mach number at a fixed pr results in a decrease of the growth rates by a factor 
of three to four up until the Mach number equals M ». For Mach numbers greater than A/., the growth rates 
of all modes level off and then with increasing Mach number, those of the slow modes begin to increase 
while those of the fast modes approach a limiting value. However, even at Mach 10, the growth rates of 
the slow modes are still small compared to those at low subsonic speeds. This, combined with the fact that 
the unstable waves have little dispersion, is a possible mechanism responsible for the observed increase in 
the flow stability. 

In this paper we report results of a study of the stability of a reacting compressible mixing layer. In 
section 2 we give the basic equations governing the mean flow and the small amplitude disturbance equa- 
tions. The boundary conditions and the numerical method are also discussed in this section. Section 3 con- 
tains a presentation of our results and conclusions are given in section 4. 

2. The Mean Flow. The nondimensional equations governing the steady two dimensional flow of a 
compressible, reacting mixing layer which lies between streams of reactants with different speeds and 
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temperatures are given by (Jackson and Hussaini, 1988) 


(P u) x + (pV) y = 0, 1 = p T, 

(2.1a,b) 

p(UU x + VU y ) = oit/,v 

(2.1c) 

p(UT x + VT y ) = (pT y ) y + (y-\)M 2 \iUy + pQ, 

(2. Id) 

p(UC it + VC iy ) = (pCjy - pA » = 1,2, 

(2. le) 

O = Dp a * M Cf C 2 ^ -e/r . 

(2. If) 


In these equations the x axis is along the direction of flow, the y axis normal to the flow, U and V are the 
velocity components in the x and y directions, respectively, p is the density, T the temperature, and C\ and 
C 2 the mass fractions. The reaction is assumed to be irreversible and of Arrhenius-type. The viscosity 4 is 
assumed to be a function of temperature. The other nondimensional parameters appearing above are: 


P = QCi^W/C P TU^i ~ AtM 

9 = E/RT„ 

D _Jo_ 

~ £/. W1 

P, = (0, - A,)W7(*1 - AiW 

M = UJa„ 


Heat release per unit mass fraction of reactant, 
Activation energy, 

Damkohler number, 

Parameter involving stoichiometry, 

Mach number, 


where a and b are the reaction orders of Cj and C 2 , respectively, the stoichiometric coefficient for species 
j appearing as a reactant, A; the stoichiometric coefficient for species j appearing as a product, 9 the sum of 
the stoichiometric coefficients of the reactants, Wj the molecular weight of species j, W the average molec- 
ular weight, E the dimensional activation energy, R the universal gas constant, B the preexponential con- 
stant in rate expression, a„ the speed of sound referred to T M Q the chemical heat release per unit mass, y 
the specific-heats ratio, and finally C p the specific heat at constant pressure. The equations were nondimen- 
sionalized by the freestream values T„, p„, (/«,, C ]>c » for the temperature, density, velocities and mass frac- 
tions, respectively. Lengths are referred to some characteristic length scale of the flow. We have 
assumed unit Prandtl and Lewis numbers in writing down these equations. The assumption of unit Lewis 
number allows us to consider linear combinations of ( 2 . Id) and ( 2 . 1 e) to eliminate the source term, which 
then admits similarity-type solutions. For a simplified Hydrogen-Oxygen reaction, typical values of the heat 
release parameter p as a function of the temperature of the moving stream are given in Table 1. Note that 
as the temperature of the moving stream is increased, P is decreased. 
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We assume that the stream at +°° is moving, while that at -» 
ditions consistent with (2.1) are 

T= U = C i — 1 , C *2 = 0 us y — > °°, (2.2a) 


C 2 _ 

T = p T , U = 0 , C x = 0 , C 2 = -r— as y - «. (2.2b) 

Ci.~ 

If is less than one, the stationary gas is relatively cold compared to the moving stream, and if p r is 
greater than one it is relatively hot 

The enthalpy for species 1 and 2 is defined by 

Hi = T+ -£-C„ i - 1,2. (2.3) 

Pi 


Combining (2.1d) and (2.1e), and using the boundary conditions (2.2), one finds 

= p r + (1 - p r + p )U + Xj-M 2 U{\-U), 


(2.4a) 


H 2 = Pj- + pr 1 + (1 - Pr - PrV + 


(2.4b) 


where § is the equivalence ratio, defined by 


P1/P2 ’ 


(2.5) 


which is the ratio of the mass fraction Cj in the moving stream to the mass fraction C 2 in the stationary 
stream divided by the ratio of their molecular weights times their stoichiometric coefficients. Note that if 
<J> = 1, then the mixture is said to be stoichiometric, and if <t> > 1 it is C, rich, while if <t> < 1 it is C i lean. 

A thin diffusion flame exists within the mixing layer and is characterized by near-equilibrium condi- 
tions; Cj = 0 on one side of the flame and C 2 = 0 on the other. In the limit of infinite activation energy this 
thin diffusion flame reduces to a flame sheet and from the relations (2.4), we see that the flow can be 
described by 

C, = 1 - (1 + <|> _1 )(1— LO* C 2 = 0, (2.6a) 


T = 1 - (1 - p T - P^Xl-t/) + ^M 2 U(1-U), 


(2.6b) 


for y>yp and 


C, = 0, C 2 = -j^MT 1 - (1 + 


(2.7a) 
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T = p r + (1 - p r + p)t/ + 3j-M 2 U(\-U), (2.7b) 

for y < y/. Here, y f gives the location of the flame sheet where both reactants vanish, and T takes the adia- 
batic flame value 7}, given by 

T f = P r + (1 - p r + P )U f + -^y-^fZ/l-ty), (2.8) 

and 


defines the flame location. 

The temperature and mass fraction profiles for the steady mean flow have been obtained in terms of 
the mean flow velocity distribution, U(y). As discussed in Part I, we assume here that 

U = j (l + tanh(i\)), (2.10) 

where 1) is the similarity variable and the Howarth-Dorodnitzyn transformation, together with the Chapman 
linear viscosity law, has been used. This profile satisfies the boundary conditions 

C/ — > 1 as q — > +«o, U —> 0 as T| — » (2.11) 

The basic mean flow considered here is thus given by (2.6), (2.7), and (2.10). 

In Figure la we show plots of T(r\) for p r = 2, M = 0, <|> = 1 and for various values of the heat release 
parameter p. If P is different from zero, the temperature always has a discontinuity in its slope at the loca- 
tion of the flame sheet. Figure lb shows plots of T for p r = 2, P = 2, Af = 0, and for various values of <|>. 
When <(> = 1, the flame sheet is located at % = 0. For <J> > 1, the mixture is C x rich and the flame location 
shifts to 0. For <p < 1, the converse is true. Finally, Figure lc shows T for various values of the Mach 
number. As can be seen from (2.6)-(2.8), increasing M increases T. It should be noted that an increase in 
P at fixed M has a qualitatively similar effect on T as an increase in M at a fixed p. 

The flow field is perturbed by introducing two dimensional wave disturbances in the velocity, pres- 
sure, temperature, density and mass fractions on either side of the flame sheet with amplitudes which are 
functions of q. In addition, the flame sheet location must also be perturbed with a wave disturbance. For 
example, the pressure perturbation is 

p = U{T\)exp[i(ax - cot)], (2.12) 

with II the amplitude. Here, for spatial stability a is complex. The real part of a is the wave number in the 
x direction, while the imaginary part of a indicates whether the disturbance is amplified, neutral, or damped 
depending on whether a, is negative, zero, or positive. The frequency c o is taken to be real. Substituting the 
expression (2.12) for the pressure perturbation and similar expressions for the other flow quantities into the 
inviscid compressible equations yields the ordinary differential equations for the perturbation amplitudes. It 
is straightforward to derive a single equation governing II, given by 
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n" - 4r~ n ' - a 2 7ir - m\v-c?]u = o, 

U—c 


which is valid on either side of the flame sheet. Here, c is the complex wave velocity 

co 


c = 


a 


(2.13) 


(2.14) 


and primes indicate differentiation with respect to the similarity variable T|. The phase speed, c pH , is given 
by to / a r . For a neutral wave the phase speed will be denoted by c N . 

The boundary conditions for n are obtained by considering the limiting form of equation (2.13) as 
T| ±oo. The solutions to (2.13) are of the form 

fl -> exp(±n ± t\), (2.15) 


where 


Cl 2 + = a 2 [l - M 2 (1-c) 2 ], Q 2 _ = ^PHPt - M 2 c 2 ]. (2.16) 


Let us define c ± to be the values of the phase speed for which Q 2 ± vanishes. Thus, 

_ i 1 _ VPr 

C+_1 M’ C ~~ A 1 ' (2,17) 

Note that c + is the phase speed of a sonic disturbance in the moving stream and c_ is the phase speed of a 
sonic disturbance in the stationary stream. At 

M = M * = 1 + (2.18) 

c± are equal. In addition to the boundary conditions at T| = ±°°, we must impose the conditions that II and 
n' are continuous across the flame sheet. 

The nature of the disturbances and the appropriate boundary conditions can now be illustrated by 
reference to Figure 2, where we plot c± versus M for a typical value of p T . In what follows we assume that 
a 2 r > a 2 -. These curves divide the c r -M plane into four regions, where c r is the real part of c. If a distur- 
bance exists with a M and c r in region 1, then Q 2 + and ft 2 _ are both positive, and the disturbance is sub- 
sonic at both boundaries. In region 3, both fi 2 + and Cl 2 - are negative and hence the disturbance is supersonic 
at both boundaries. In region 2, Cl 2 + is positive and Cl 2 - is negative, and the disturbance is subsonic at +°° 
and supersonic at and we classify it as a fast mode. Finally, in region 4, Cl 2 + is negative and Cl 2 - is 
positive so the disturbance is supersonic at +«> and subsonic at -«>, and we classify it as a slow mode. 

One can now see that the appropiate boundary condition for either damped or outgoing waves in the 
moving and stationary streams are, respectively, 


n ^ 

if c r > c + . 


if c r < c + , 

(2.19a) 

n -> e a -". 

if c r < c_, 


if c r > c_. 

(2.19b) 



( 2 . 20 ) 
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To solve the disturbance equation (2.13), we first transform it to a Riccati equation by setting 


G = 


TV 

a7Tl ' 


Thus, (2.13) becomes 

OJV t* 

G' + arc 2 - [•—— - 4t]G = a[T - M z (U-c) 2 ]. (2.21) 

U-c T 

The boundary conditions can be found from (2.19) and (2.20), with G continuous across the flame sheet. 

The stability problem is thus to solve equation (2.21) for a given real frequency to and Mach number 
M, with U and T defined by (2.10) and (2. 6b), (2.7b). In addition to M, the important parameters of this 
problem are pj~, (3, and q>. The eigenvalue is the wavenumber a. Because this equation has a singularity at 
U = c N , we shall integrate it along the complex contour (-6,-1) to (r^O) and (6,-1) to (q^O), where gives 
the flame sheet location defined by (2.9). Using a Runge-Kutta scheme with variable step size, we choose 
an initial a and compute the boundary conditions from (2.19). We then iterate on a, using Muller’s 
method, until the boundary conditions are satisfied and the jump in G at Cty-,0) is less than KT 6 . All calcu- 
lations were done in 64 bit precision. 

3. Results. In all of our calculations we have taken y = 1.4; p r = 0.5, 1, 2; 0 < p S 5; $ = 0.5, 1, 2; 
and 0 2 M < 10. 

Lees and Lin (1946) have proven that if a neutral mode is to exist in region 1, the phase speed will be 
given by c N = £/(tj c ), where q,. is found from the regularity condition 

S(q) = -£(r 2 -^) = 0. (3.1) 

a q «q 

The corresponding neutral wave number, a N , must be determined numerically. The eigenfunction is called 
a subsonic neutral mode. Because f is discontinuous at ty for nonzero P, S(q) will also be discontinuous at 
this point. Fiqure 3 is a plot of S versus q for various values of p. As one can see from examination of 
this figure, S can have a single root, two roots one of which corresponds to q positive and the other nega- 
tive, or two roots one of which is a one-sided zero. The roots of S, which corresponds to phase speeds that 
are subsonic at both boundaries, are the phase speeds of subsonic neutral modes. The one-sided zero of S 
may or may not yield a phase speed of a neutral mode. Finally, there can also be singular neutral modes, 
ones whose phase speeds are not given by (3.1) and are subsonic at the boundaries. If the phase speed 
corresponding to a zero of S is supersonic at either or both boundaries it may or may not be that of a neu- 
tral mode. This can only be determined numerically. 

3.1. M = 0. The introduction of a chemical reaction, in the form of a flame sheet, has complex 
effects on the flow stability. In order to show these effects clearly, we first consider the case of zero Mach 
number and examine the variation of the eigenvalue, a, with the heat release parameter p and equivalence 
ratio <f>. In Figure 4a we show a plot of the growth rate (-a,) versus the real part of the wavenumber (a r ) 
for the case of p r = 0.5, 4> = 1, and various values of p. Note that all of the curves pass through the origin, 
which corresponds to co = 0, and also cross the a r axis, corresponding to co > 0. These give the two neutral 
waves. We classify these as true stability curves since the unstable waves are found between two neutral 
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waves. 

The curves shown in Figure 4a are for values of P <, 0.98. For this particular set of parameter values 
we were not able to find solutions to the stability problem whose eigenvalues had the behavior shown in this 
figure for any values of P > 0.983. In order to understand the reason for this we carried out an extensive 
search for eigenvalues in the a plane while varying p. We discovered that the eigenvalue problem has two 
unstable solutions for any P > 0 and the other values of the parameters. Some of these results are shown in 
Figure 4b. It is clear that there is a saddle point for P about 0.984 at a = 0.24-0.60*. For P < 0.984 there 
are two sets of eigenvalue curves. The lower set appears to be the physically relevant one because it yields 
the eigenvalue relation for the nonreacting mixing layer as P approaches zero. We believe that the upper 
branch, while corresponding to a correct mathematical solution of the eigenvalue problem, has no physical 
relevance because none of these upper branch curves intersects the a r axis and hence there are no neutral 
solutions (a, = 0) or stable solutions (a, > 0). Similarly, we believe that the solutions on the other two 
branches, the dashed branches in Figure 4b, are also non-physical because each of them crosses the a r axis 
at only one point in the right half plane, that is with a r > 0 and hence c r > 0. In other words, these solution 
branches have only one neutral solution with c r > 0 and we therefore believe that these are not physically 
meaningful solutions. Thus for this set of parameter values, the physically acceptable modes only exist if 
0 < P < 0.983. 

We have found similar behavior in the eigenvalue spectrum for other values of the parameters, that is, 
there is a critical value of P (P c ) beyond which one of the physically acceptable modes no longer exists. 
This value of P c depends on all of the other parameters of the problem, and the general trends are similar. 
As Pj- is increased, p c also increases. Typical values of p c are given in Table 2. The implications, and pos- 
sible causes, of the existence of p c will be discussed in section 4. 

In Figures 5, 6, and 7 we show the variation of the neutral phase speeds c N with p for various values 
of Pr and <j>. First, recall that for the nonreactive case (Part I), with p = 0 and <J> = 1, only a single neutral 
mode exists at M = 0, given by (3.1), and is classified as fast or slow depending on whether c N > 1/2 or 
c N < 1/2, respectively. We will show that the addition of heat release (P > 0) can cause both fast and slow 
modes to coexist, although both types may not always exist for all values of p. Also, there may exist singu- 
lar neutral modes, those which are found numerically but do not satisfy the regularity condition (3.1). This 
is in direct contrast to the results of Part I. 

It can be shown from (3.1), and from (2.6b), (2.7b) and (2.10), that fast waves only exist for 

P S 1 - Pt4> 

with corresponding neutral phase speed 

_ P + P^ 

N P + (1 + PrM ’ 

while slow waves only exist for 


(3.2) 


(3.3) 


P > P^ - 1 


(3.4) 
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with corresponding neutral phase speed 


Pr 

Cn ~ l + p + p r ‘ 


(3.5) 


Note that the value of the slow neutral phase speed is independent of the equivalence ratio 4>. Also note 
that there is a mode switch located at P = 0 and 4> = p r _1 ; i.e., the modes interchange their characteristics at 

P = 0- 


The phase speeds of the neutral modes for p 7 = 0.5 are plotted in Figure 5. Note that there are both 
fast and slow subsonic neutral modes whose phase speeds are given by (3.3) and (3.5). In addition, there 
are both fast and slow singular subsonic neutral modes, adjancent to their corresponding regular neutral 
modes, whose phase speeds are independent of p. It is clear from the results shown here that it is the slow 
modes which cease to exist when P > P c and that p c increases with increasing values of the equivalence 
ratio <(>. The fast modes, on the other hand, only exist for those values of p which satisfy (3.2). From Fig- 
ure 5c one can see that there will be a mode switch for ()> > 2 at P = 0. The phase speeds of the neutral 
modes for Pr = 1 are plotted in Figure 6. Again we see that there are both fast and slow regular neutral 
modes with adjacent singular modes. Note that there is a mode switch at p = 0 as <p is increased past one. 
That is, for <)> < 1, only the slow modes exist at P = 0, and conversely for $ > 1. The critical value p e , 
beyond which there are no acceptable solutions, only affects the slower branch. Finally, Figure 7 shows the 
variation of the phase speed with P for Pr = 2. The mode switch occurs when exceeds 0.5. It is clear 
that P c is greater than 5 for all values of <|> for which results are shown here. 

In addition to the neutral modes there are, of course, associated unstable modes. The maximum 
growth rates of the unstable modes for various values of Pr and <t> are plotted as a function of p in Figures 
8-10. It is important to note the changes of scale between these figures. 

In Figure 8 the variation of -a^ with P is shown for p 7 = 0.5 and <> = 0.5, 1, 2. The slow modes 
are much more unstable than the fast modes for those ranges of P for which the slow modes exist, i.e. 
P < p c . The maximum growth rates of these slow modes first decrease slightly and then increase with an 
increase in p. The maximum growth rates of fast modes increase with <(> for any p but are much smaller 
than those of the slow modes. Similar results are shown in Figure 9 for p T = 1. The switching between fast 
and slow modes past $ = 1 which was noted in our discussion of the phase speeds of the neutral modes is 
apparent in this figure. The growth rates of the fast waves increase with <(>, while those of the slow waves 
decrease. Finally the plots of the maximum growth rate versus P for p r = 2, shown in Figure 10, have the 
same general behavior as the two previous figures. That is, there is a mode switch at p = 0 past <|> = 1/2. 

The results shown in these figures can be summarized as follows. An increase in Pt- causes a decrease 
in the maximum growth rate of the most unstable modes. At fixed p 7 and $ the growth rate of the slow 
modes increases as P becomes large, while the growth rate of the fast modes approaches a limiting value. 
As we shall see, this is the same generic behavior as results from increasing the Mach number while hold- 
ing all other parameters fixed. Thus, for sufficiently large p, the slow modes are the most unstable ones. 
For fixed p, an increase in (j> causes an increase in the maximum growth rates of the fast modes, while 
decreasing the maximum growth rates of the slow modes. Finally, as in Part I we have found that the phase 
speeds c ph of the unstable modes lie between the phase speed of the regular neutral mode and its adjacent 
singular neutral mode. 
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3.2. M > 0. In this section we show the variation of the phase speeds of the neutral waves and the 
maximum growth rates of the unstable waves as a function of the Mach number for various combinations of 
P T , P, and <)>. 

In Figure 11a we show the phase speed of the neutral modes for the nonreactive (P = 0, <|> = 1) mixing 
layer (Part I) for p r = 1/2. From this figure, taken from Part I, one can see that there is only a single sub- 
sonic neutral mode in region 1. This mode crosses over the sonic curve at M s , the Mach number at which 
the phase speed equals that of a sonic wave, and is transformed into a supersonic neutral mode in region 4. 
In addition, a fast supersonic neutral mode appears at M* in region 2. In regions 2 and 4 there are unstable 
modes with phase speeds between that of the supersonic neutral mode and that of the sonic neutral mode. 
Thus, there are two bands of unstable frequencies for Mach numbers greater than M*. The band in region 2 
is a group of fast and the band in region 4 is a group of slow unstable modes. The phase speeds of both the 
fast and slow modes have a small range about the average, so that little dispersion of wave packets is 
expected, with a reduction in the dispersion as the Mach number is increased. 

It can be seen from Figure lib that chemical reaction has a major effect on c N . First, recall from Fig- 
ure 5b that the slow subsonic neutral mode only exists for p < p c at M = 0, while the fast subsonic neutral 
mode only exists for P > 1/2. As can be expected, the value of P c is also a function of the Mach number. 
Thus, we find that the slow neutral mode only exists if the Mach number is greater than the critical value 
M c (see Figure lib). Apart from this the values of c N on this branch are only slightly affected by changes 
in P and the effect is quite small at large M. Finally, one can see that a fast subsonic neutral mode is 
present for all M because P is greater than 0.5. Note that the values of c N for this fast mode in region 1 
increase markedly with increasing p. This is also true for the lower Mach number range in region 2, but all 
of the c N curves for different P appear to asymptote to a single curve for large M. Associated with the fast 
and slow modes there are singular subsonic neutral modes in region 1. The range of phase speeds of the 
unstable modes in region 1 lies between the fast modes and their associated singular mode, and between the 
slow modes and their associated singular mode. In regions 2 and 4, the range of phase speeds of the 
unstable modes lies between the fast (slow) supersonic neutral modes and the associated fast (slow) sonic 
neutral modes. The range of the phase speeds of the unstable modes increases with 3 and is quite large for 

Mach numbers around M», thus yielding an increase in dispersion. However for Mach numbers much larger 
than M. the dispersion is essentially independent of p. 

Figure 12 shows the phase speed of the neutral modes for p T = 1, <J> = 1 and various values of P as a 
function of the Mach number. The results for p = 0 are shown in Figure 12a and those for P > 0 in Figure 
12b. The curves of c N versus M are symmetric about the line c N = 1/2 because the mean velocity profile is 
symmetric about T| = 0 (where U = 1/2) and so is the temperature profile for any p if <J> = 1. Hence, the 
subsonic neutral mode in region 1 of Figure 12a splits into a symmetric pair of fast and slow supersonic 
modes at M*. From Figure 12b at can be seen that this is also true for any nonzero p. The only exception 
is for sufficiently large p, where the slow neutral modes only exist if the Mach number exceeds M c . For 
this case, the singular subsonic mode in region 1 has a phase speed of 1/2, and the range of phase speeds of 
the unstable modes lies between this mode and the corresponding fast and slow neutral modes. As P 
increases the phase speeds of the fast (slow) neutral modes increase (decrease) at any fixed M. Again it 
should be noted that all of the curves for both the fast and slow neutral modes at different P are asymptotic 
to a single curve for large M. As for the previous case, an increase in P causes an increase in the range of 
the phase speeds of the unstable modes, in the vicinity of M ., and hence an increase in the dispersion. 
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Finally, in Figure 13 we show similar plots of c N versus M for p r = 2. In Figure 13a, with p = 0, 
there is only a subsonic fast neutral mode in region 1 which is transformed into a supersonic neutral mode 
in region 2. The slow supersonic neutral mode only exists in region 4. The effects of increasing p are 
similiar to those of the cases discussed above. There is an increase in the phase speed of the fast modes as 
P is increased. The phase speed of the slow supersonic modes of region 4 decreases with increasing p and, 
for p > 1, a singular subsonic neutral mode with c N = 1/2 appears in region 1. Correspondingly, for P > 0, 
there also exists a singular subsonic neutral mode in region 1 associated with the fast modes. Again the 
curves for the fast and slow neutral modes are each asymptotic to a single curve for large M. As in the pre- 
vious two cases, an increase in P causes an increase in the range of the phase speeds of the unstable waves 
and hence an increase in the dispersion. 

We have also carried out calculations of c N versus M for fixed p T and P with various values of 4>. 
The results can be easily summarized. As <t> increases the phase speed of the slow mode is unchanged, con- 
sistent with the Mach zero results (see equation 3.5). The only effect is a change in the critical value of M c 
below which this neutral mode does not exist; the smaller the value of <|> the larger is the value of M c . The 
phase speed of the fast neutral modes decreases with increasing <>. Finally, all of these curves appear to be 
asymptotic to a single curve as M -» °°. 

In addition to the phase speeds we have calculated the maximum growth rates of the unstable modes 
for these values of p r , P, and <|> as a function of the Mach number. The results are shown in Figures 14-16. 
Again, it is important to note the changes of scales between these figures. 

The maximum growth rate as a function of M with p r = 1/2, <|> = 1 and various P is plotted in Figure 
14a. For P = 0 an increase in the Mach number from zero results in a decrease in the growth rate of the 
slow mode by a factor of three to four up to A/., and for higher Mach numbers the growth rate levels off 
and eventually begins to increase with increasing Mach number. An increase in P causes the growth rate of 
the slow modes to increase, but the variation of the growth rate with M is similar to that for p = 0. How- 
ever it must be recalled that, if P > 0, the slow modes only exist for Mach numbers greater than M c , and 
this is reflected in the curves for p > 0.983. The growth rates of the fast modes are much smaller and 
increasing p has only a slight effect on them. Figure 14b contains similar results, but with P = 2 and vary- 
ing <|>. An increase in <J> causes an increase in the maximum growth rate of the slow modes and a decrease 
for the fast modes at low and moderate Mach numbers. At higher Mach numbers, the growth rates are 
essentially independent of <|>. 

Similar results, but with p r = 1, are given in Figure 15. It can be seen that the effects of varying p 
on the growth rate (Figure 15a) is somewhat more complex than for the previous case. The maximum 
growth rate of the slow modes at low Mach numbers first decreases and then increases with increasing p, 
consistent with the Mach zero results of Figure 9. At higher Mach numbers, increasing P causes a mono- 
tonic increase in the maximum growth rates of the slow modes. Again there is a three to four fold decrease 
in the growth rates of the slow modes as the Mach number increases from zero to M,, and a slight increase 
at higher Mach numbers. Apart from the appearance of the fast modes at low Mach number, increasing p 
has little effect on their maximum growth rates. Finally, in Figure 15b, an increase in 4> with fixed p results 
in an increase in the maximum growth rate of the fast modes and a decrease in the maximum growth rate of 
the slow modes at low and moderate Mach numbers. At higher Mach numbers, the growth rates are again 
independent of <J>. 


- 13 - 


In Figure 16 similar results are given for p r = 2. These results are consistent with the Mach zero 
results of Figure 10. In particular, Figure 16a shows that the slow modes are strictly increasing and those of 
the fast modes are strictly decreasing with increasing p at zero Mach number. Thus, increasing p causes a 
substantial decrease in the maximum growth rates of the fast modes and a substantial increase in the max- 
imum growth rates of the slow modes at low and moderate Mach numbers. As in the previous case, the 
results of Figure 16b show that increasing <)> yields an increased in the growth rates of the fast modes while 
decreasing those of the slow modes. As before, the growth rates of both are independent of <J» at higher 
Mach number. 

The results presented in Figures 14-16 thus show that the maximum growth rates of the slow modes 
increase as the Mach number becomes large, while that of the fast modes approaches a limiting value. As 
we have seen, this was the same generic behavior as results from increasing the heat release parameter p at 
zero Mach number. 

Figures 17 thru 20 are plots of selected two dimensional neutral eigenfunctions for p r = 1, <J) = 1, 
Mach 5 and increasing values of p. These plots show the variation of n with T] r on the contour t|, = -1. 
All of these have been normalized so that the maximum of the absolute value of n is unity. The eigenfunc- 
tions shown are all fast supersonic neutral modes. Note the rapid variation of both the amplitude and phase 
near = 0. Because these modes are fast, they show exponential decay in the subsonic region and oscilla- 
tions with constant amplitude and linear phase in the supersonic region. As P is increased, the variation of 
the amplitude near T) r = 0 increases markedly. Also, the decay of the amplitude in the subsonic region is 
more rapid and the rate of change of the phase in the supersonic region increases. The behavior of the 
eigenfunction with increasing heat release, P, is similar to the behavior with increasing Mach number. 

4. Conclusions. The addition of combustion in the form of a flame sheet has important, and complex, 
effects on the flow stability. In contrast to the nonreacting case, we have shown the existence of multiple 
regular and singular subsonic neutral modes in region 1. For Mach numbers greater than A/., there are two 
bands of unstable frequencies (just as in the nonreacting case); one a group of fast supersonic modes and the 
other a group of slow supersonic modes. In general, an increase in p causes an increase (decrease) in the 
phase speed of the fast (slow) neutral modes. Since the range of phase speeds is increased at low and 
moderate Mach numbers as the heat release parameter p is increased, there is an increase in the dispersion 
of wave packets. Finally, we have found that the addition of chemical heating has almost no effect on the 
phase speeds at higher Mach numbers. 

The maximum growth rates of the unstable modes decrease by a factor of three to four as the Mach 
number approaches M,, even with the presence of the reaction. Just as in the nonreacting case, for Mach 
numbers greater than Af., the growth rates level off and those of the slow modes eventually begin to 
increase with increasing Mach number, while the growth rates of the fast modes approaches a limiting 
value. The same behavior results from increasing P while the Mach number is held fixed. Finally, if the 
stationary gas is colder than the moving gas, p r < 1, the effect of increasing the equivalence ratio <t> is to 
increase the growth rate of the slow modes and decrease the growth rate of the fast modes. On the other 
hand, if p r > 1, increasing <)> has the opposite effect. However, at higher Mach numbers, changes in <}> have 
little effect on the growth rates for any p r . Finally, an examination of the eigenfunctions shows that 
increasing P has the same generic effect as increasing the Mach number. 
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As was discussed in section 3.1, we found that with p r = 1/2, 4> = 1, and M = 0 there was a saddle 
point in the complex a plane for p about 0.984. For any p > 0 there were two branches of a. As was dis- 
cussed in that section, we accepted only that branch which yielded the eigenvalue relation for the nonreact- 
ing mixing layer as p approached zero. Similar behavior was also found for other values of p r and Mach 
numbers, but generally the saddle point appeared at larger values of p. Thus, one can conclude that there 
will be a value of P beyond which one of the branches ceases to have physically meaningful solutions to the 
stability problem. We suspect that this is due to the use of the flame sheet approximation. In this model, 
although T is continuous everywhere, f is discontinuous across the flame sheet and the magnitude of the 
jump in T increases with p. If we are correct, the activation energy 0 must be taken to be large but finite. 

This study is the only comprehensive study, of which we know, of the inviscid spatial stability of a 
reacting compressible mixing layer. We do not know how sensitive our results are to the assumptions used 
in this study. In particular, we have assumed unit Prandtl and Lewis numbers, used Chapman’s linear rela- 
tion between viscosity and temperature, and approximated the mean velocity profile by a hyperbolic tangent. 
In addition, we have taken the limit of infinite activation energy which reduces the diffusion flame to a 
flame sheet Despite these limitations, we believe that this systematic study is an important first step in 
classifying and understanding the complex effects that chemistry has on the stability of compressible free 
shear layers. As mentioned above, the next step is to consider a more realistic model of the chemistry and 
the thermodynamics. This will then yield mean velocity, temperature, and mass fraction distributions which 
will be continuous and have continuous derivatives across the flame. However the calculation of the mean 
field as well as the perturbation solution will be more difficult since now the velocity, temperature, and 
mass fraction equations are coupled. We have begun this study with large, but finite 0, and hope to report 
the results at a later date and will compare those results to the benchmark results reported here. 
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C P 

Qm/^ l w l 

P 

(K°) 

(cal/mole K°) 

(kcal/mole) 


1000 

9.03 

516.5 

3.2 

1500 

9.94 

516.5 

1.9 

2000 

10.60 

516.5 

1.4 


Table 1 . Typical values of the heat release parameter (3 as a function of 
temperature. 


<b 

Pr 

0.5 

1.0 

2.0 

0.5 

0.7 

0.984 

1.8 

1.0 

2.0 

3.1 

>5 

2.0 

>5 

>5 

>5 


Table 2. Typical values of p c at M - 0 as a function of |3 r and <j>. 


















Figure 11a. Plot of neutral phase speeds (solid) and sonic speeds (dashed) versus Mach number for 
(Jr = 0.5. 0 = 0, and = 1. 









Figure 14a. Plot of maximum growth rates of the fast and slow modes versus Mach number for 
Pr=0J, p = 0, 1, 2, 5, and + = 1. 
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Figure 15a. Plot of maximum growth rates of the fast and slow modes versus Mach number for 
pr=l, P = 0, 1,2,5, and * = 1. 



Figure 15b. Plot of maximum growth rates of the fast and slow modes versus Mach number for 
p T = 1, P = 2, and + = 0.5, 1, 2. 
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Figure 17. Plot of the two dimensional fast supersonic neutral eigenfunction n(t|) along the contour 
r) = t| r - i. The solid curve corresponds to the amplitude and the dashed curve to the phase. M = 5, 
pr = 1, $ = 1, P = 0, with (% = 0.184813, a* = 0.215661, c N = 0.85696. 


1.0 




Figure 18. Plot of the two dimensional fast supersonic neutral eigenfunction n(T|) along the contour 
- i. The solid curve corresponds to the amplitude and the dashed curve to the phase. M = 5, 
p T = 1, <(> = 1, p = 1, with to* = 0.235413, a* = 0.269485, c N = 0.873565. 




Figure 19. Plot of the two dimensional fast supersonic neutral eigenfunction ITCn) along the contour 
n = Tl,- i. The solid curve corresponds to the amplitude and the dashed curve to the phase. M = 5, 
3r=l, <►=!, 3 = 2, with cd* = 0.276698, a* = 0.312299, c N = 0.886001. 



Figure 20. Plot of the two dimensional fast supersonic neutral eigenfunction I1(T|) along the contour 
ti = T\ r - i. The solid curve corresponds to the amplitude and the dashed curve to the phase. M - 5, 
3r = 1, <|> = 1, 3 = 5, with co„ = 0.356124, a N = 0.390535, c N = 0.911887. 
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